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ABSTRACT 



For the purpose of studying the effect of varying quality upon the 
pressure drop accompanying isothermal, two phase, two component flow, 
a horizontal, one-inch flow channel of circular cross section, into which 
air could be injected at various points, was designed and constructed. 

Preliminary tests were made using water and air injected at a single 
position. The resultant data confirmed the findings of Chisholm and 
Laird for turbulent flow in a smooth tube. 

Tests were then conducted in which the quality, the percent mass of 
gas phase of the total mass flowing, was varied by injecting known amounts 
of air at seven injection points along the 42- foot test section. The 
Chisholm- Laird correlation provided pressure drop predictions to within 
25% of the measured values. 
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lo Introduction. 



Accurate prediction of the pressure drop experienced by two phases 
of a substance (or substances) flowing concurrently through a pipe or 
channel is required in many engineering situations. Some common uses of 
two phase flow are: pipeline gathering systems in the gas and oil industry; 

heat exchangers; horizontal tube evaporators; steam generation lines; re- 
actor coolant channels; pipeline contactors in chemical plants; and conden- 
sate lines in refrigeration plants. 

The term, ' two phase flow", used herein, will mean either two compon- 
ent, two phase flow, or single component, two phase flow, according to the 
context . 

As fluids pass through horizontal, closed channels or pipelines, they 
experience a loss of pressure through energy dissipation due to turbulence 
and friction between the walls and the fluid. For single phase flow the 
pressure drop may be predicted with adequate precision using the Darcy- 
Weisbach equation and a friction factor obtained either from the ioody 
Diagram or, if the pipe is smooth and the flow turbulent, calculated from 
the Blasius formula. Unfortunately, a simple summation of the pressure drops 
due to each of the phases flowing separately does not predict the actual loss 
of pressure. 

The increased complexity of computations encountered is due to one or 
more of the following reasons. 

The hydraulic diameter of each of the phases is variable and generally 
undefinable as flow progresses. The cross section available for flow of 
one phase is reduced as the second phase is introduced. The unsteady, in- 
tensely turbulent flow, which is characteristic of high mass flow rates in 
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two phase flow, is an efficient energy absorber. Much energy is irrecover- 
ably expended creating ripple, wave, plug, and slug flow. No longer is 
the flow boundary a relatively smooth pipe wall, having become Instead a 
rough, mobile, irregular Interface between the two phases. In view of such 
problems, considerable effort has been required to obtain good correlations. 
The history of two phase flow investigation extends back to 1939 when 
Boelter and Kepner /18/* first attempted to obtain a workable correlation. 

In 1942 Dittus and Hildebrand /19/ presented a method which used both mathe- 
matical equations and graphical solutions for the determination of pressure 
drop in oil vapor-oil mixtures flowing through furnace coils. Then in 1944 
came the first generally accepted correlation for the prediction of two 
phase flow pressure drop, that of Martinelli /3/. During succeeding years, 
Martinelli and his associates attacked the problem of pressure drop during 
the forced circulation boiling of water /5/, and produced a revised cor- 
relation in which the liquid fraction (the fraction of the cross-sectional 
area of the pipe occupied by liquid) was integrated into the prediction /4/, 
In 1954 Alves /20/ reported experimental results within 15% of the pre- 
dicted pressure drop, and classified the seven types of flow that were 
observed. White and Huntington /7/ in 1955 proposed a new correlation 
which brought the predicted values to within 12% of the measured values, 
for the very restricted type of flow known as "ripple flow". In 1955 
Chenoweth and Martin /6/ presented a correlation that came within 20% for 
turbulent flow, Chisholm and Laird /!/ modified Martinelli'' s flow- rate- 
ratio parameter, introduced the effect of roughness, and obtained a cor- 
relation to within 15% without recourse to a graphical solution, A two 
phase friction factor developed by Bertuzzi, Tek and Poettmann in 1957 /21/ 

/ / denotes reference in Bibliography. 
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provided the basis for a correlation which produced a standard deviation 
of 20.87o. In more recent years. Hatch and Jacobs /22/ investigated the 
prediction of pressure drop in the two phase flow of hydrogens, and Kordy = 
ban /23/ presented a flow model that provided a logical explanation for 
the applicability of the Chenoweth-Martin correlation /6/, 

A preliminary library search produced no reports which dealt with the 
pressure drop associated with flow of varying quality. Consequently a 
horizontal, one-inch, flow channel of circular cross section, into which air 
could be injected at various points, was designed and constructed. 

Preliminary tests were made using water and air injected at a single 
position. The resultant data confirmed the findings of Chisholm and Laird 
for turbulent flow in a smooth pipe. 

Tests were then conducted in which the quality, the percent mass of 
gas phase of the total mass flowing, was varied by injecting known amounts 
of air at seven injection points along the 42 “foot test section. The 
Chisholm* Laird correlation provided pressure drop predictions within 25% 
of the measured values. 
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2. Description of Equipment, 



a. Flow Channel Apparatus . The 42-foot test section^ Figs. 1 and 2 ^ 
were made of ext ruded, acrylic plastic tubing of a nominal one-inch inside 
diameter. Seven sections of the tubing, six feet in length, were joined 
in station connector blocks (Fig. 3) using 0-ring seals. This construc- 
tion provided very satisfactory sealing with a minimum reduction in the 
smoothness of the tube. The actual inside diameter varied from 0,995 to 
1,025 inches both in the diametral plane and along the tube. 

An effective diameter of 1.013 inches was used. By assuming the 
tube to be smooth, using equations (1) through (8) (listed in Section 3-b) 
with the values obtained from the single phase water runs, and treating the 
diameter, D, as the unknown variable, the resultant equation was solved for 
the effective diameter. This diameter was then used, as required, in sub- 
sequent calculations. 

Figure 4 displays the results obtained from the single phase water 
runs using this diameter. It is evident from Fig. 5 that the single phase 
air runs suffered no adverse effects from the use of the diameter thus cal- 
culated. Included in the effective diameter were the effects of the eight, 
station connector blocks. Slight discontinuities in the smoothness of the 
tube were introduced by the two tubing ends butting against an insert of the 
same tubing glued inside the block, and by the annular series of holes for 
each of the pressure taps. Distances between stations are shown in Appendix 
I. 

b. Auxiliary Equipment. An overall view of the test equipment is pro- 
vided in Fig. 1, A more detailed view of the auxiliary equipment shown in 
Fig. 6, In this photograph the following pieces of equipment are evident: 
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Number 



Item 



1. Water Drain Stop Cock 

2. Rotameter Selector Switch Valve 

3. Rotameter Selector Switch Block 

4. Injection Air Rotameter 

5. Water Rotameter 

6. Differential Manometer Inlet Pulsation Damper 

7. Station Pressure Manometer Pulsation Damper 

8. Manometer Line Bubble Trap 

9. Water Main Inlet to Bubble Trap (with air bleeder valve 
directly above) 

10. Precision Potentiometer 

11. Electric Timer 

12 Timer Control Box 

20. Cathetometer 

21. Manometer Bank for Measuring Rotameter Pressure 

22. Differential Mercury Manometer 

23. Air Pressure Regulators in Series 

24. Differential Water Manometer 

Where these pieces of equipment are indicated on the g(S:neral apparatus 
schematics Fig. 7* the numbers correspond. 

The main features of the auxiliary equipment^ the selector switcheSg 
may be seen in Figs, 8 and 9. In these figures most details of their 
construction are apparent. Two details not evident are^ 

1. Whenever gluing of the plastic was required^ a mixture 
of thin acrylic shaving in ethelyene dichloride was inviscid enough to 
penetrate the thinnest crack or abutment of the pieces being joined. This 
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method provided joints of more than adequate strengthe 

2o The aluminum alloy (24ST) slide valve^ item 2j, Figo 9 
had slots milled along half of its length^3/16 inch steel tubing laid 
thereinj, and was returned to its cylindrical shape by filling with 
cement. This eliminated the difficult task of successfully drilling small 
diameter holes for a long distance in a relatively soft material. The holes 
which are evident in the valve j lead directly to the tubing connections at 
the right hand end of the valve. These In turn are connected to the inlet 
and outlet of a rotameter. 

The rotameter selector switch, Fig. 9, allowed insertion of the ro- 
tameter into each station’s air supply line successively, only momentarily 
disrupting the continuity of any of the supply lines, and the "’steady^* state 
of the flow. 

The station pressure selector switch, Fig. 8, was used to transfer the 
manometer pressure leads from one station to the next, and could also be 
used for easily obtaining differential pressures between successive pairs 

of stations. 

The bubble trap, item 8, Fig. 10, kept air bubbles from entering the 
manometer pressure lines. Water was admitted to the top of the bubble trap 
through the valve, item 9. Entrapped air was then bled off the top of the 
trap as it accumulated. Station pressures were measured in all cases by 
mercury or water manometers vented to the atmosphere. All differential 
pressures between stations were measured with a mercury or water differen- 
tial manometer. In the case of differential pressures, all liquid levels 
were measured to within 0.05 millimeters using a Gaertner cathetometer . 

Temperatures were measured using copper-constantan thermocouples refer- 
red to an ice bath. The wiring diagram for this circuit is indicated in 
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Fig. 11. Temperature measurements made at the inlet (thermocouple #1) 
and at the outlet (thermocouple #2) indicated that the flow was isother- 
mal to within plus or minus 0.1®F. 

The high pressure air source fluctuated between 80 and 100 psig. 

This fluctuation was reduced to plus or minus 0,1 inches of water by the 
use of two pressure regulators connected In series. 

The water main pressure was found to be sufficiently constant, varying 
in mass flow rate by only 0.007 Ibm/sec, when the tests were conducted 
after hours of normal water usage. 

The weigh tank scale was checked with dead weights. 

Rotameter readings were converted to cubic feet per minute using the 
curves supplied by the manufacturer. 

The test section was horizontally aligned, using a K&E Tilting Level, 
to within plus or minus 0.03 inches, A chalk line, snapped between the 
extremities of the test section support, provided a straight reference 
line for longitudinal alignment. See Fig, 2, 

Figs. 12 and 13, with the wiring schematic. Fig. 14, show the method 
employed to measure the time required for a set mass of water to flow. 
Movements of the tank switch contact arm caused this mass of water to vary 
from 148 to 154 Ibm, and consequently the "set” mass frequently required 
redetermination. 
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3. Procedure. 



a. General . Except as modified in subsequent sectionsj, each test 
run was conducted in a similar manner. Temperatures were stabilized by 
running both the water and air for an hour before the test. Atmospheric 
temperature and pressure were recorded. The quick closing valve on the 
weigh tank was closed and the weighing proceeded automatically. Water and 
air injection temperatures were recorded. The overall pressure drop^, station 
pressures^ and air rotameter pressure were then noted. 

b. Single Phase Qualifying Tests . In order to establish the validity 
of measurements made, tests were conducted to determine the friction pres-^ 
sure drop for single phase water flow and single phase air flow. The re- 
sults of these tests are displayed in Figs. 4 and 5. In all cases in 

this investigation, data reduction was accomplished using a digital com- 
puter. 

The determination of friction factors, Reynolds numbers, and pressure 
drops for the case of a fluid flowing alone was made using the following 
equations : 



( 1 ) 



= yf K. j 



where : 




U/l = Mass flow rate, Ibm/sec, 

3 

= Density of liquid, Ibm/ft 




= Average velocity of flow, ft/sec, the velocity 
of the liquid which would occur if only the 
single phase were flowings 



rearrangement gives, 
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Using this velocity, the Reynolds number was calculated by: 



(3) , 

where : - Reynolds number, dimensionless, 

D = Effective diameter of tube, ft, 

= Absolute viscosity of liquid, Ibm/ft-sec. 
Since the test section was a smooth tube, the Blasius formula provides 
the applicable friction factor: 

_ .3Ua 



(4) 



/ = 



/■ 



Friction factor, dimensionless. 



A/r*® 

One form of Bernoulli’s equation states: 

<=> ^ ■ 23,^ 

In the case of horizontal flow with constant area and fluid density, ^ 
equation (5) reduces to the form: 









(6) 



i^-1 i 









^ - ~ir ■ 

where : Af 5s Pressure drop due to liquid flowing alone in 

the tube, . ' 

ss Friction head, ft-lbf/lbro. 



The Darcy-Weisbach equation states: 




L = Length 


Equating (6) and (7): 




or. 


je ■ ^ 


(8) = 


APu = ^ 


L 






■H 



where : AR s Pressure drop per unit length, psf/ft. 

c. Two Phase Flow . 

1. Background of the Chisholm and Laird Correlation . 

Use is made in this paper of a two phase flow pressure drop correlation 
developed for smooth and rough tubes by D. Chisholm and A. D. K. Laird 
/!/ in 1956. The reasons for choosing this correlation are twofold; the 
data could be correlated to within 20%, and the portion of the correla- 
tion used, was completely in mathematical terms, which facilitated data 
reduction by a digital computer. 

When correlations which were available to Chisholm and Laird were 
used, the percent error in prediction of pressure drop increased as a 
given set of flow conditions deviated farther from the conditions of flow 
in a smooth tube. These errors were commonly of the magnitude of 30%. 

To produce a better correlation Chisholm and Laird collected data 
from concurrent, horizontal, isothermal flow of water and air through 
smooth and rough tubes of approximately one- inch bore. The ranges of ex- 
perimental data which they correlated, and the types of tubes used, are 
summarized, with those of the present investigator, in Tables 1 and 2 on 
the following pages. 
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TABLE 1 



RANGES OF DATA 



QUANTITY 


CHISHOLM & LAIRD 


1 AUTHOR 


UNITS 


Mass velocity 
water 


39 to 600 


51 to 236 


2 

Ibm/ sec-ft« 


air 


0.1 to 20 


0.589 to 2.10 


Ibm/ sec- ft 


Reynolds number 


water 


4500 to 80,000 


5790 to 26,000 




air 


1000 to 140,000 


4090 to 14,580 




Temperature 


water 


59 to 74 


52.3 to 61.8 


degrees F. 


air 




64.8 to 68.3 


degrees F. 


Viscosity 


-5 - 5 


-5 -5 




water 


77x10 ^to 63x10 - 


85x10 to 73x10 ^ 
1.2x10“ to 1.21x10' 


Ibm/ft-sec 
Ibm/ ft-sec 


air 


1.2x10“ to 1.23x10“ 



Tube Surfaces 

pipe roughness ratio 

Acrylic plastic* 0.000 

Smooth brass** 0.000 

Commercial galvanized tube** 0.0025 

Brass tube, concrete Internal surface** 0.013 

Brass tube, Internal thread** 0.028 (measured) 

0.037 (apparent) 

Brass tube, non-unlformly distributed sand**0.045 
Brass tube, uniformly distributed sand** 0.068 

General Remarks 

Flow In both phases was considered turbulent throughout the ranges 
considered. Mean arithmetic values of the water temperature at the Inlet 
and the outlet were used when determining water properties. Air properties 
were determined with no appreciable error by assuming that the air was dry 
and the temperature quickly became the same as the mean water temperature 
upon mixing. 



* Author* s tube 
** Chisholm & Laird’s tube 
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TABLE 2 



SPECIFIC RANGES OF MASS FLOW RATE 
FOR EACH TYPE OF FLOW 
STUDIED IN THIS INVESTIGATION 

AIR, Ibm/sec 

0,00330 to 0,0118 

0.00348 to 0.00950 

0.000426 to 0,00142 
(per injection point) 



TYPE OF FLOW WATER, Ibm/sec 

Single phase water 0,508 to 1.32 

Single phase air - - - - 

Two phase single Injection 0.286 to 1.27 

Two phase multiple 0.694 to 1.30 

injection 
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Through Introduction of the effects of roughness into the Martinelli 
correlation Chisholm and Laird arrived at the following correlation of 
their data. They define a flow-rate-ratio parameters X s 



where; Liquid-mass velocity based on tube cross sections Ibm/sec-ft 

2 

= Gas-mass velocity based on tube cross sections Ibm/sec-ft g 
Absolute viscosity of liquids Ibm/sec-ft, 

=s Absolute viscosity of gas, Ibm/sec-ft, 

= Density of liquid, Ibm/ft^, 

3 

jO = Density of gas, Ibm/ ft . 

The pressure-drop parameter, * is then correlated by equation 

arl 

APl X 

where ; Total (friction) pressure drop over an increment of length 

Friction pressure drop for liquid flowing alone in tube, psf/ft, 

X 



^ ^ Vl^PLP I 



vn= A 
A« 



= Constants for a particular liquid flow rate and tube surface. 



Reynolds number where the liquid flows alone, 
\ = Pipe friction factor for rough tube, 

X§ =» Pipe friction factor for smooth tube. 

Equation (10) may now be expressed as 



( 11 ) 









^7 
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The liquid fraction, Rl , is correlated with the pressure-drop para- 
meter, ^ as follows; 

AFL 



( 12 ) 



A/tp 



0.6 



(for the smooth tube). 



where: ^ = Liquid fraction, the area occupied by the liquid 

divided by the cross-sectional area of the tube* 

A correlation of the liquid fraction, Ru , with the flow- rate-ratio para- 
meter, ^ , for the smooth tube is: 



(13) ^ I — ■ t — . (for the smooth tube), 

X X 

li. Single Injection Tests . After the single phase tests were 
completed, a series of runs was made in which water and air were Injected 
upstream of the test section. The correlated pressure drop was obtained 
in the following manner; 

1. X determined using equation (9). 

2. Ru was obtained from equation (13). 



3. f::' J-i could then be calculated from equation (12), 

AR. 

The results of these single injection tests, together with the results 
of the multiple injection tests (described in the next section), are pre- 
sented graphically in Figs, 15, 16, and 17, where the coordinates used are 
those found to be the most useful by Chisholm and Laird, 

iii. Multiple Injection Tests , Having thus shown that the ap- 
paratus is capable of producing data comparable to that presented by 



14 



Chisholm and Laird, test runs were conducted in which air was injected 
in known amounts at seven points along the test section. 

The purpose of this multiple injection of air was to provide portions 
of the test section with two phase flow whose liquid fraction was able to 
be determined. In addition, it was desired that the liquid fraction vary 
from injection point to injection point. The design of the test section 
allowed injection of air at each station connector block. Since the amount 
of air injected at each station was measured, the liquid fraction was calcu- 
able, and constant for each portion of the test section between injection 
points. The correlation process, set forth previously in section 3-c-ii, 
was then used and the pressure drop found for each constant- liquid- fraction 
section. The pressure drops were then summed over the whole test section, 
and the overall pressure drop obtained. The results of these multiple in- 
jection tests are displayed in Figs, 15, 16, and 17. 
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4. Results and Conclusions, 



A horizontal channel suitable for multiple air Injection Into the 
primary flow of water was designed, constructed, and tested. 

The results of single phase tests, Figs, 4 and 5, Indicated that 
the test section behaved as a smooth tube. Consequently, the friction 
factor computed by the Blaslus formula was the only one employed during 
subsequent reduction of the two phase runs. 

These Figures show that the pressure drop accompanying two phase flow 
of varying quality may be predicted to the same degree of accuracy as 
constant quality two phase flow, by a suitable selection of the "constant 
quality" sections. A logical next step would be the prediction of pressure 
drop In a two phase, varying quality, Isothermal situation, such as a 
saturated steam generating tube, by applying the correlation to sections 
of the tube small enough to be of "constant quality", from one end to the 
other. This could be done on a digital computer and the results compared 
with experimental data. A further extension of this line of Investigation 
might well start with mounting the variable quality test section vertically. 
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APPENDIX I 



LOCATION OF PRESSURE TAPS ALONG THE TEST SECTION 



Station 

Number 


Distance to next 
station^in Inches 


Distance to station indicated 
from station Onegin inches 


1 


71.7 


0.0 


2 


71.72 


71.7 


3 


71.75 


143.42 


4 


71.8 


215.17 


5 


71.75 


286.97 


6 


71.8 


358.72 


7 


71.7 


430.52 


8 


0.0 


502.22 



Distance of station 1 from inlet --37.22 inches 
Distance of station 8 from outlet --35.53 inches 
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APPENDIX II 

DRAWINGS, PHOTOGRAPHS, AND' GRAPHS 



21 



II 

r 




22 



Fig. 1 General View of Equipment 
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Fig. Z Overall View of Test Section 
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Construction Details of a Station Connector Block 
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REYNOLD'S NUMBER 




Fig. 4 Friction Factor vs. Reynolds Number 
Single Phase Water Flow, 
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REYNOLD'S NUMBER 




Fig. 5 Friction Factor vs, Reynolds Number 
Single Phase Air Flow, 
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Fig. 6. Auxilliary Equipment Detail 
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FIGURE 7 PIPING ARRANGEMENT OF TEST APPARATUS 
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Rotameter Station Selector Switch Detail 
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Fig. 10. Control Panel Detail 
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Fig. 12. Tank Switch Detail 




Fig. 13. Timer Control Box 
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FIGURE 14. 
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Fig. 15 Experimental Pressure Drop Ratio vs. 

Reciprocal Liquid Fraction in Turbulent 
Flow of Air and Water. 
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•Fig. 16 Experimental Pressure Drop Ratio Minus One 
vs* Chisholm-Laird Flow Parameter in Tur- 
bulent Flow of Air and Water. 
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Pig, 17 Experimental Pressure Drop *Ratlo vs. Predicted, 
Pressure Drop Ratio in Turbulent Flow of Air 
and Water. 
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